基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

矩的显式积分算法研究与应用

矩的求解通常被用于求解有限元、体积、惯性矩等问题中.基于矩的叠加性,首先给出了在三维空间中计算域的离散方式,并推导了矩的显式积分公式,随后将其推广到n维空间中,该表达式易于在计算机上实现;设计了矩的并行计算算法,并通过Fortran和Python混编的方式,实现了矩的并行计算;对多重精度下的样例数据给出了一个算例,实现了零阶矩和二阶矩的计算,并和串行算法、逐次降维算法作出比较,进行了效率分析和误差分析.结果显示,矩的显式积分并行计算算法易于程序实现,并且在效率上高于串行算法,能够很容易推广到高维空间,该算法具有高度可并行性,误差主要来自计算域离散.      

悬停状态倾转旋翼/机翼干扰流场及气动力的CFD计算

基于一套高效通用的多层运动嵌套网格技术,建立了适合悬停状态倾转旋翼机旋翼/机翼气动干扰特性分析的高效混合计算流体力学（CFD）方法。在倾转旋翼/机翼贴体网格区采用可压雷诺平均Navier-Stokes（RANS）方程作为主控方程,过渡/背景网格区选用Euler方程,湍流模型选用Spalart-Allmaras模型。时间推进上采用双时间推进格式进行非定常求解,并在方法中运用了SPMD（Single Program Multiple Data）模式的并行加速技术。在此基础上,首先分别采用UH-60A直升机旋翼及XV-15倾转旋翼机旋翼作为数值算例,验证了CFD方法的有效性。然后着重对倾转旋翼/机翼的非定常干扰流场及气动力分布特征进行了数值研究,模拟得到与实际情况相符的“喷泉效应”干扰现象。计算结果表明,干扰作用使得倾转旋翼相对于孤立旋翼拉力数值减小了3%,但总的拉力系数损失达到了17%,证明悬停状态下气动干扰对飞行器气动性能有重要影响。     

航天器着陆碰撞动力学CUDA并行计算技术

航天器着陆接触碰撞过程中的瞬态响应计算,涉及到几何非线性、接触非线性和材料非线性,计算过程复杂,计算量大。为此,文章将统一计算架构(CUDA)技术应用到接触碰撞问题的并行显式积分计算,研究了质量矩阵和刚度矩阵的并行生成方法,并将基于中心差分法的求解过程并行化实施。仿真算例表明:并行算法计算结果准确,较大地提高了计算仿真分析的效率,可为航天器着陆过程中冲击动力学计算的快速和高精度实现奠定基础。     

过渡状态倾转旋翼气动力模拟的高效CFD方法

为显著减少倾转旋翼过渡飞行时气动力CFD模拟的计算代价,提出并建立了适合倾转旋翼过渡状态气动特性分析的高效混合CFD方法。首先,提出了适合于过渡状态模拟的嵌套网格系统,并发展了相应的挖洞和贡献单元搜寻方法。在此基础上,结合叶素理论和动量理论建立了旋翼气动力模拟的简化虚拟桨盘模拟方法(Virtual rotor model,VBM)。为了能够分析倾转旋翼气动力的细节特性,多层运动嵌套网格系统和单指令多数据流(Single program multiple data,SPMD)并行技术被引入来建立精确的旋翼模拟方法(Real blade model,RBM)。然后,通过将VBM和RBM方法结合,构建了适合倾转旋翼过渡状态气动特性分析的高效(Hybrid blade model,HBM)方法。最后,通过对有试验值对比的悬停状态典型旋翼和7A旋翼分别验证了VBM和RBM方法的有效性。分别采用3种方法预测了过渡状态不同倾转角下旋翼的气动特性,VBM表现出最优的计算效率,能用于倾转旋翼总体气动性能的分析。HBM方法在保证流场求解精度的基础上,相对于高精度的RBM方法节省了1/3的计算时间。     

C/E复合材料螺旋铣孔切屑形状与切削温度研究

碳纤维增强环氧树脂基体(C/E)复合材料具有优异的性能,广泛用于制造飞行器结构件。在加工C/E复合材料连接孔过程中容易出现分层、撕裂等加工缺陷,而切削温度过高是导致加工缺陷的主要原因之一。螺旋铣孔作为新的制孔技术,加工C/E复合材料时表现出许多优于传统钻孔的特点,受到广泛关注。本文采用专用加工装备开展了C/E复合材料制孔试验,利用红外热像仪实时监测切削温度,分析了螺旋铣孔加工参数对切削温度的影响规律。根据螺旋铣孔运动学规律和切削原理,从未变形切屑形状特点入手分析了螺旋铣孔切削温度的来源和影响因素。研究结果可对深入理解螺旋铣孔加工机理、制定合理的螺旋铣孔加工工艺提供依据。     

多导线对称分布原子磁分束器和干涉仪设计分析

提出了采用多根对称分布载流导线构成原子分束器的方法,包括三导线和四导线磁导引。阐述了原子分束器的分束机制,用Ansoft Maxwell 2D软件计算给出了部分载流导线移动到不同位置时的磁场分布图。通过分析了导引中心的变化,发现只要通过改变载流导线之间的相对位置就可以来实现从单路到双路导引和三路导引的转换。由此分别设计由3根对称分布载流导线构成的原子双路分束器和4根对称分布载流导线构成三路分束器,然后运用Monte Carlo方法模拟验证其原子分束功能。最后以三线对称分布分束器为基础构成Mach-Zehnder原子磁干涉仪。     

全机模型静强度分析的多区域并行计算

在Windows XP 64位操作系统环境下,利用若干普通PC机搭建了机群系统,实现了分布内存式并行（DMP）计算。采用3节点和4节点薄壳单元和2节点梁单元建立了约260万自由度（DOF）的全机有限元模型。采用封装的METIS区域分解方法（DDM）,并利用大型通用有限元软件MSC.Marc,实现了全机静强度模型的并行求解,为大型复杂有限元模型的求解提供了一种高性价比的解决方案。     

一种三维结构网格并行算法

基于三维结构网格,以MPICH2作为并行开发环境,研究了对等模式下流场的并行计算问题。在三维结构网格已经生成的前提下,针对计算区域的划分,网格输出格式和数据通信,提出了一种并行算法。算法的核心思想是将需要通讯的其他区域的单元按区号从1到n排列,并获得本区域通信单元的I-J-K编号与其他区域通讯单元的排列序号的映射关系。最后,在IBM BladeCenter HS22刀片机平台上进行了算例验证,结果表明算法具有很好的可靠性与通用性。